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Abstract 

In this article, we address the dynamics and bifurcations of a wide class of stochastic differential 
equations around singular points where both the drift and diffusion functions vanish. We apply 
these results to stochastic versions of the pitchfork, Hopf and saddle-node bifurcations. Ac- 
cording to the Holder coefficient of the diffusion function around the singular point, we identify 
different regimes, and the multiplicative noise case appears as a singular transition. 
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1. Introduction 

A large class of phenomena can be described using non-linear ordinary differential equations 
(ODEs), in such distinct domains as economics, physics or biology. Many analytical methods 
have been developed to characterize the solutions of such systems. In particular, the widely ap- 
plied theory of bifurcations allows a qualitative characterization of the asymptotic regimes of 
nonlinear dynamical systems. This theory ensures that the number or stability of fixed points and 
periodic orbit can change according to a few generic scenarii, and the system at the bifurcation 
point can be locally mapped to a simple polynomial dynamical system, called normal form [1,2]. 
This theory is very handy to study nonlinear differential equations and is almost systematically 
used to characterize the asymptotic behavior of complicated systems. Investigating the impact 
of noisy perturbations on such systems is hence of great interest, and is currently an active field 
of research. In the last decades, it has been shown in many different areas that applying noise 
on a non-linear system can lead to many counter-intuitive phenomena, such as noise-induced 
stabilization (3), oscillations or transport (5] |6) and stochastic resonance ||7]|8). From a 
mathematical perspective, understanding the interplays between non-linearities and noisy pertur- 
bations is a great challenge, with many applications. Several tools have been introduced, ranging 
from the theory of random dynamical systems (RDS) J9][l0), the study of moment equations 
ifTTI fT2ll or multiscale stochastic methods for slow-fast systems jT3l[T4l . Unfortunately no sys- 
tematic method in the flavor of bifurcation theory for the analysis of the dynamics of nonlinear 
ODEs exist, and this is the central problematic of the present manuscript. 
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We focus here on the dynamics of stochastic differential equations (SDE) in which both the 
drift and the diffusion functions vanish at some particular point, called singular point. This type 
of equations may appear for example in diffusion approximations of population models IT5l or 
in the modeling of noisy neural conductances [16]. The question we address is how the interplay 
between the drift and the shape of the diffusion function affects the behavior of the system. An 
important contribution of RDS theory to the field of stochastic bifurcations is to distinguish be- 
tween phenomenological (P) bifurcations, characterized by a change in the shape of the invariant 
measure, and dynamical (D) bifurcations, characterized by a sign change of the stochastic Lya- 
punov exponent. We investigate the questions of stochastic stability and bifurcations (P and D) in 
a comprehensive manner, using combinations of appropriate tools such as Lyapunov functionals 
and exponents, stationary solutions of Fokker-Planck equations and Feller scale functions (see 
e.g. IfTTij for an extensive account on these concepts). Mathematically, the choice of a vanishing 
diffusion coefficient provides a framework to study subtle competitions between the drift and the 
noise, leading to a rich and generic phenomenology. 

With the aim of characterizing the behavior of the solutions of SDEs with drift given by 
the normal form of usual bifurcations, we will start by analyzing a general one dimensional 
stochastic differential equation 

dx t = f(x t ) dt + y(x t ) dW, 

where W = (W t ) t>0 is a Brownian motion defined on a complete probability space (Q, T , P) 
endowed with the natural filtration {^F,) t of W. Notice that several classical normal forms can be 
reformulated in this setting and our general results will be applied to the cases of the pitchfork, 
Hopf or saddle-node bifurcations in section [3] We assume that there exists a singular point x* 
such that fix*) = y{x*) = 0. The system has a trivial solution distributed as a Dirac measure 
localized at the singular point x*, and we shall address the stochastic stability of these solutions. 
In contrast to the case of ODEs, there are at least three different notions of stochastic stability: 

1. The almost sure exponential stability defined by the property 

limsup - log \x t — x* | < a.s. 

/— >oo t 

for any initial condition Xq e R. The almost sure exponentially instability defined by the 
property 

liminf - log \x, - x*\ > a.s. 

for any initial condition Xq e R (see e.g. ifTSl ). This definition implies a strong notion 
of stability close to the usual exponential stability property of deterministic dynamical 
systems: almost all trajectories exponentially fast converge to or diverge from x*. In our 
stochastic context, this notion can appear quite restrictive, and we will hence also discuss 
a weaker notion of stability: 

2. The stochastic stability (or stability in probability), defined by the property that for all 
p e (0,1) and r > there exists 6 depending on p and r such that 

¥({\x,-x*\<r Vt>0})>l-p 

whenever \xq — x*\ < 6 and otherwise it is said to be unstable in probability. The solution 
x* is stochastically asymptotically stable (or asymptotically stable in probability) if it is 
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stable in probability and for every /u G (0, 1) there exists 770 depending on p such that 

P(lim|x,-JC*| = 0) > 1 -fi 

whenever |jcq - x*\ < 770. 
3. The stability (resp. instability) in the first approximation of a stationary solution X II 1 71 
Chap. 7.1], notion widely used in the RDS theory (9), consists in that the identically 
zero solution of the linearized equation around X is almost surely exponentially stable 
(resp. unstable). This criterion will be chiefly used when dealing with non-Dirac stationary 
solutions. 

To be more specific, we now consider that x* = and assume that f(x) is a smooth function, 
which is often the case in applications and a basis of the deterministic bifurcation theory. In 
detail, we assume that f(x) — A x + g(x) where g(x) = o(x) at x = (i.e. /'(0) = A), and hence 
consider the following class of SDEs on R: 

dx t = (Ax t + g(xt)) dt + y(x,)dW t (1) 

with y(x) ~ cr\x\ a at x = 0. The behavior of the solutions will be described as a function of the 
local behavior of g and y around similarly to the deterministic case, and these properties will 
precisely characterize the stability of the solutions. However, in contrast with the deterministic 
case, the existence of stationary distributions will depend on the properties of g and y at infinity. 
Intuitively, this is due to the fact that stationary solutions can reach values far from the singular 
point, where the local behavior of the solution is no more relevant. In Figure[T[ the main question 
of the paper is schematically described, in the setting of a particle in a double well potential 
subjected to noisy fluctuations that vanish on the saddle point and increase as a function of 
the distance. According to the shape of the noise function y (mainly the parameter a), the 
particle will have two distinct fates : either it falls into the well and keeps fluctuating without 
ever reaching the saddle point, or it eventually reaches it after some finite time and remains steady 
for all time. 

The paper is organized as follows. We start by analyzing in detail in section [2] the behavior 
of the solution of a general diffusion equation ([1]). This section contains several abstract results 
summarized in Theorem[T]and a general approach that will be useful for the analysis of classical 
bifurcations. We investigate the existence, uniqueness and blow up of the solutions and show 
that noise can prevent explosion. The stochastic stability of the singular point is characterized 
and the absorption time at the singular point is shown to be almost surely finite if and only if 
a < 1. In that case, we show that quasi-stationary distributions exist when a < |. In the 
case a > 1, we study the existence of other stationary distributions. Beyond the mathematical 
developments, this paper is also resolutely intended for engineers, physicists and other scientists 
interested in the qualitative behavior of noisy systems described through stochastic differential 
equations. Therefore, readers interested in applications may skip this section and jump directly to 
section[3]where they will find readily applicable results. In that section we build on the results of 
the general analysis to analyze the stochastic pitchfork, Hopf and saddle-node bifurcations with 
power diffusion function. Interestingly, in all the cases, any system locally reducible around an 
equilibrium point to one of these normal forms displays the same qualitative behavior, provided 
that the drift satisfies an additional condition at infinity. This property is analogous to the well- 
known universal unfolding property for deterministic dynamical systems. In all three cases, we 
put in good use the characterization of the dynamics by developing applications to the Wilson 
and Cowan system lfl9ll . a simple neural population model widely used in neurosciences. 
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noise vanishes at 




Figure 1 : Stability of a singular saddle point in double-well potential perturbed by power diffusion coefficient. A particle 
located at the saddle point remains steady. If a small perturbation is applied, two scenarii are possible according to the 
local behavior of y(x) (see theorem^: either noise brings back the particle to the saddle point (left) or the particle falls 
into the well and fluctuates without ever reaching the saddle point (right). 



2. General Theory 

In this section we analyze the behavior of the solutions of a general equation ([TJ under a few 
assumptions on the function g and y: 

(HI). g(x) ~ jjx l+K at x — 0, for some /i e R and k > 

(H2). g(x) < -vx l+,} when \x\ > A for some A > 0, for v e R (generally assumed non-negative) 
and p > 0. 

(H3). y(x) ~ cr\x\ a at x — 0, for some cr > and a > 

(H4). y(x) ~ d x s at x = +oo with d > and 5 > 0. 

(H5). y(x) only vanishes at x - 

(H6). Either (i) v > or (ii) v < and 6 > 1 + /3/2. 



Assumptions (HI) and (H3) describe the behavior of the solutions close from zero and will gov- 
ern the stability of the singular point 0. The assumptions |(H2)||(H4)| and |(H5)| are global proper- 
ties that will be chiefly used for explosion matters, and existence of stationary distributions. Such 
assumptions are generally not necessary in the deterministic case: indeed, in the local unfolding 
of a bifurcation, the presence of a non-zero fixed point implies the local existence of a manifold 
of fixed points under some non-degeneracy conditions. However, in the stochastic case, as soon 
as the singular point (where noise vanishes) is unstable, the system will visit very regions very 
remote from the singular point with positive probability, and these latter conditions will come 
into play. 

The stochastic process distributed as a Dirac at zero (X, = a.s.) is solution of the equation 
whatever the parameters. Studying the stationary solutions and their stability will make essential 
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use of the differential operator _£ associated with Ito's representation of the diffusion acting on 
twice differentiable functions V e C 2 (R, R): 



£V(x) = (Ax + g(x)) V'(x) + -y(x) 2 V"(x). 



Beyond the well-posedness of the equation studied below in section 2.1 the main results of 
the present section can be summarized in the following: 



Theorem 1. Under the assumptions \( HI )\ \( H5)\ we have: 

• for j < oc < 1, the singular point is stable in probability and any solution is absorbed 
infinite time almost surely at the singular point. If moreover a < 3/4, there exists an 
exponentially attractive quasi-stationary distribution for the trajectories that did not reach 
the singular point. 

• for a — 1, the singular point is stable for A < o~ 2 /2 and unstable otherwise, and two 
stationary solutions exist for A > cr 2 /2. 

• for a > 1, the singular point is stable for A < and unstable otherwise, and two stationary 
solutions exist for A > 0. 



2.1. Existence and uniqueness of solutions 

The first question arising with equation ([1} is its well-posedness. Indeed, this equation does 
not satisfy the standard conditions for existence and uniqueness of solutions: the drift function 
is neither globally Lipschitz-continuous, nor holds a linear growth condition. Moreover, the dif- 
fusion coefficient is also singular: for a < 1, the diffusion coefficient is not Lipschitz continuous 
at 0, and for a > 1 it does not satisfies the linear growth condition generally imposed for the 
existence and uniqueness. Based on existing results [20 1, we prove the following property: 

Proposition 1. For any a e [j* 00 )* there exists a unique strong solution to the equation ([TJ. 



Under the assumptions (H2) (H4) and (H6) the solution never blows up infinite time. 



Proof. In the case where a > 1 it is clear that both the drift and diffusion functions are locally 
Lipschitz continuous. This property implies that the strong uniqueness property of the solutions 
holds (see e.g. EOl Chap. IX.2]). Moreover, the drift and the diffusion functions are locally 
bounded, ensuring the existence of solutions up to an explosion time. These two properties 
ensure strong existence and uniqueness of solutions up to an explosion time. 

The case a < 1 is slightly more complex. In that case, the linear growth condition is clearly 
satisfied, but the diffusion coefficient is not Lipschitz at 0. However, the method presented in 1201 
theorem IX. 3. 5] based on the properties of the local time at zero of the diffusions ensure that we 
have the existence and uniqueness of the solutions of ([T| for a > j and non-uniqueness of 
solutions for a < \ . 

Let us now address the non-explosion property. To this purpose, we use Feller's test for 
explosion (see EH)- The function Ax+ y s ^i +l is locally integrable on both the intervals (-°o, 0) and 
(0, oo). We consider the interval (0, oo) for instance (the same principle applies for the interval 
(-oo, 0)), and for some c > 0, we define the scale function: 



The process almost surely never diverges before possibly reaching zero as soon as soon as p(x) 
is diverges at oo. At infinity, we have: 

for some s > 0, arbitrarily small when x is sufficiently large. This quantity is lower bounded at 
infinity in the case where v > or in the case 6 > 1 + /?/2, implying in both case that p(x) — > oo 
when x — > oo. We hence conclude that the process does not blow up in finite time in these cases, 
which ends the proof. □ 

Remark. This property illustrates a first stabilization effect of power diffusions, namely the blow-up 
prevention. In the case where the solutions to the noiseless equation do not blow up in finite time, power 
diffusion coefficients cannot create a stochastic blow up. In the case where the solutions to the noiseless 
system blow up in finite, a power diffusion function prevents blow up provided that the diffusion coefficient 
diverges fast enough at infinity compared to the diffusion function. 

2.2. Stability of 5q 

We now address the stability of the equilibrium 6q, and prove the following: 
Theorem 2. Let us denote by 6q the stationary solution X, — a.s. for all t. We have: 

• For a < 1, <5o is asymptotically stable in probability whatever /ieM and cr 0. 

• For a — \, 6q is: 

(i) . asymptotically stable in probability for A < Moreover, if g(x) — -jj.x 1+K for 

[M,k > and y(x) — crx, the stationary solution <5o is almost surely exponentially 
stable for X<^. 

2 

(ii) . unstable in probability for A > ^- 

• For a > 1, <5o is : 

— stable in probability for A < 

- unstable in probability for A > 



Proof. Case a < 1: Let us consider V(x) — x a as a Lyapunov functional for the dynamics. This 
function is clearly C?(H), i.e. it is twice continuously differentiable except at x = 0. Moreover, 
we have: 

£V = ax"' 1 (Ax + g(x)) - a( l ~ a) 7 (x) 2 x a - 2 

Since we assume that a < 1 and x = o(g(x)) at 0, the leading term close to is of order 
-a(l - a) /2 cr 2 x 3a - 2 which is strictly negative for sufficiently small x. More precisely, there 
exists r > such that V e C^iU, ) where U r is the open ball of radius r and such that £V < for 
all x e U r . Moreover, for any e such that < s < r we have V(r) > and £fV < -c E < for 
all r > x > e. Theorem V.4. 1 of [ 17 1 therefore applies and concludes the proof of the stability in 
probability of the solution do whatever the parameters A e R, cr > and a G [1/2, 1). 
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Case or = 1: We start considering the case A < a /2 and define rj — ^ - -4 > 0. The 
Lyapunov function V : jc h-» jc ,; satisfies: 

XV (jc) = 7?x"(i + - l)y 2 (x)) + 77X"- 1 



Under assumption (HI) we have jc^(jc) ~ /j.x 2+i< , hence is equivalent, close from zero, to 

/ (>7-l) y 2 (x)\ 1 cr 2 

_j = -(A- T )^. 

Since (/l-cr 2 /2) < 0, the same argument as in the case a < 1 ensures asymptotic stability of 6o in 
probability for A < cr 2 /2. If we now further assume g(x) = -fix l+K and y(x) — crx, considering 
V(x) = jc 2 , we have: 

• £y (x) = (2A + cr 2 )x 2 -2 /J x 2+K < (2A + cr 2 )V(x), 

• \V'(x)cr x\ 2 = \2ctx 2 \ 2 = 4o- 2 V(x) 

Using the stochastic stability theorem 1 18, Theorem 3.3], we deduce that: 

1 Act 2 - 2(2A + cr 2 ) cr 2 
limsup - log \X(f)\ < - A — — almost surely. 

t— >oo t 4 2, 

By definition of the almost sure exponential stability, we conclude that the solution is almost 
surely exponentially stable for any A < 

Let us now deal with the case A > cr 2 /2. We define V(x) = -log(jc). This function tends to 
infinity as x — > and belongs to C^R). Moreover, 

2x x 

which is equivalent at zero to -A + cr 2 /2 < 0. Hence there exists r > such that £V < for 
any x < r, and moreover we have for any s > and r > x > e, V(x) > and £.V(x) < —c s < 0. 
Theorem V.4.2 of IfTTl applies and ensures that the solution 6q is not stable in probability. 
Case a > 1: For V(x) = x 2 around zero, we have: 

£V = 2Ax 2 + 2xg(x) + y(x) 2 

and under the assumption that g(x) ~ fvc 1+K at zero and a > 1, the term £y is equivalent to 2Ax 2 
and hence locally has the sign of A. For A < 0, we can directly apply theorem V.4.1 of IfTTl . For 
A > 0, we use V(x) = - log(x) again and conclude that is not stable in probability. □ 



2.3. Extinction time 

We know that X, — a.s. is a stationary solution of the SDE. It is an absorbing state of the 
diffusion, i.e. as soon as the solution reaches the singular point, the solution is almost surely iden- 
tically stuck at the singular point for the whole subsequent evolution. We are now interested in 
the extinction probability, namely the probability of reaching zero in finite time. We distinguish 
between the cases a < 1 and a > 1 . 
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2.3.1. Case a < 1 

Proposition 2. Under assumptions 



WTWH6) 



the first hitting time of zero of the solution of 



dx t = (Ax + g(x)) dt + y(x)dW, 
is almost surely finite for any a < 1 and cr > 0. This is also the case for v — and A < 0. 

Proof. In order to establish this property, we use a sufficient condition derived from Feller's test 
for expositions and evidence the almost sure finite exit time of the interval ]0, °°[. The finiteness 
of the exit time of ]0, oo[ will readily imply the finiteness of the hitting time of 0, since the process 
is almost surely bounded for all finite time as proved in proposition [TJ under the assumptions of 
the proposition. 



Let us define the interval / = (0, oo), and let c e / arbitrarily chosen. Under assumption (H5) 
we have that for all x e I, the diffusion coefficient is always strictly positive, hence the function 
(1 + \Ax+g(x)\)/y(x) 2 is locally integrable on /. The finiteness of the exit time of the process from 
the interval / can be established through the analysis of Feller's scale functions. For convenience, 
we introduce the function: 

7 2 (0 

Feller's scale functions written as a function of G read: 



P(x) d = £ exp {-2G(y)} dy 
m{dx) ± f 2t *$<™> dx 

V W = f*(p( x )-p(y))tn(dy) 

We prove that (a) lim^oo p(x) = oo and (b) lim M o + v(x) < oo. The property (a) has been 
addressed in the proof of proposition [TJ We hence only need to prove the property (b) which is 
slightly more subtle. The function v(x) is indeed defined as an integral, which is well behaved 
except possibly at x = 0. Close to this point, we have: 



G(y)-G(z)= " , T' dt 



A -(y 2(l - a) -z 2{l - a) )+-^-(/-z^) 



2o~ 2 (l — a) cr 2 $ 

with (p = 2 + k - 2a and this function is smooth and bounded close to zero. Let us denote by M 
the supremum of \G(y) - G{z)\ on [0, c], and by F{x,y) the function defined for < y < x < c: 

F(x,y)= f exp(2(G(y)-G(z)))dz<e 2M (x-y). 

Jy 



It is then clear that: 

Jc r (y) Jc 7 00 

which is integrable at zero since the integrand is of order x l ~ 2a and hence integrable at zero under 
the condition that a < 1. By application of Feller's test for explosion ||2T1 Prop. 5.32 (iii)] we 
conclude that if S denotes the exit time of the process from the interval /, we have W(S < oo) = 1. 
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Moreover, we have seen in proposition [T] that under the assumptions of the proposition, the 
process is almost surely bounded. We hence conclude that necessarily the process reaches in 
finite time and is absorbed (this is also a direct consequence of Feller's test for explosion ||2T1 
Prop. 5.22] under the condition that p(x) — > oo when x — > oo if p(x) > -oo when x — > + , which 
is a simple consequence of the above estimates). 

Let us now deal with the case v = and A < 0, we have G(x) - /bc 2(I ~ a) + C for some constant 
C. The same reasoning as above ensures that p(x) tends to infinity as x — > oo. Demonstrating 
that v(x) is bounded at zero did not use the fact that v + and hence applies in the present case. 
We conclude that P(S < oo) = 1, hence the process almost surely reaches in finite time. 

□ 

We hence proved that for a < 1, the solutions almost surely hit zero in finite time, where they 
stay trapped for the subsequent dynamics, as soon as the diffusion coefficient cr is nonzero. 

When cr = 0, this is never the case: in the case A < 0, the solutions of the related ODE 
converge exponentially fast towards zero but never actually reach it, and in the case A > 0, the 
solutions escape from zero exponentially. It is hence of interest to analyze the behavior of the 
system in the small noise limit. We focus on the case A > (in which case is unstable in 
the noiseless system and a.s. reached in finite time in the stochastic case) and use Freidlin and 
Wentzell large deviations theory lt22l to analyze in detail the behavior of the system in the small 
noise limit. In order to obtain precise quantitative estimates, we will assume that g(x) = —p.x l+K 
with k > and fi > and y(x) = cr\x\ a , corresponding to a double-well potential. Note that the 
method developed here does not depend on the choice of g(x). 

We analyze the properties of the extinction time of the system with initial condition chosen 
as the fixed point of the noiseless system: 



r :=hvNf>0; X, = | X 



l/K 



In the case a < 1, one can expect from the large deviation theory that IE [to] becomes exponen- 
tially large as cr — > 0, i.e behaves in the leading order as e^°~ where V is a positive constant. 
More precisely, we have: 



Proposition 3. For all A < a < 1 and A > 0: 



lim cr 2 log IE [to] = A 

cr— >0 



2(1-ci)/k 



1 1 



2(1 -a) 2(1 -<*)+* 



> 



Proof. By application E2l Theorem 4.1], the problem is reduced to compute at the boundary 
x = the quasipotential U(x), solution of the Hamilton-Jacobi equation: 

(Ax - fix l+K )U'(x) + —(U'(x)) 2 = (2) 

with boundary condition !/((-)«) = 0. That is, one has to compute: 



U(x)= £ du (3) 
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Simple algebra yields the expression: 



1/(0) = A 



2(\-a)/K 



1 



1 



2(1 -a) 2{\-a)+K 



(4) 



□ 



As a consequence, for a < 1 fixed, the extinction time becomes exponentially large, which 
explains the transition between the regime cr > (where the system a.s. hits zero) and cr — 0. A 
critical boundary of this domain is the point or = 1: in that case, the system a.s. never hits zero 
in finite time as shown in section l2~.3.2l 

In particular, proposition[3]ensures that the convergence as a — > 1~ with A > fixed displays 
a super-exponential time behavior: 

lim lim cr 2 log IE [to] = oo 

ml- cr-tO 

Going deeper into the analysis of the above convergence allows to understand more precisely 
what happens in the critical parameter region where both a — > 1~ and A — > + , in the limit of 
small noise cr — > 0. Such a "blow up" reveals the following picture around the point (a, A, cr) = 
(1 , 0, 0) in parameter space: 

Corollary 1. Assume that A — > + and a — * 1~ such that 

A 



(I -a) 

Then, one has the following three regimes: 

• if c — 0, the extinction time is sub-exponential : 



c 6 [0,oo] 



(5) 



lim lim cr log IE [r ] = 

a-»M-»0er-»0 



//0 < c < oo, the extinction time is exponential : 



lim lim cr log E [to] = — ^ — rr 

a ->U-,0 ->0 5 UJ 2p 2( - l - a ^ K 



if c — oo, the extinction time is super-exponential : 



lim lim cr log IE [to] = oo 

a-tU->0o--t0 



where lim is constrained by Eq. (Bl). 
<r-»M->0 jj 

Proof. This proposition is a consequence of the analysis of the different asymptotic regimes of 
the critical exponent c* obtained in proposition [3] 



2(1-o0/k 



1 



1 



2(1 -a) 2(1- a) +* 
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The leading term in the considered limit is proportional to A l+2( - 1 a ^ K /(l -a). A series expansion 
provides the following expression: 

c * ~ V^(^TT + l Al °^ + OWog(Af(l - a)). 

The different regimes identified in the corollary immediately follow from the analysis of the 
limits of this exponent as A — ♦ + and a — > 1~ subject to the condition (BJ. □ 

Remark. Note that the regimes do not depend on p. or k qualitatively and the property is purely local: as A 
goes to zero, the integration range becomes increasingly small. The property is hence valid for equation {T} 
with function g(x) equivalent at zero to -px [+K . However, the asymptotic exponential divergence rate quan- 
titatively depends on the local behavior of the flow close to zero, as observed in the case < c < oo. 

From a dynamical systems viewpoint, this type of solutions corresponds to a relatively sharp 
convergence towards the equilibrium. Indeed, the dynamics stops at the almost sure hitting time 
of the singular point and the whole subsequent dynamics is deterministic and constant. This will 
not be the case for more regular diffusion coefficients (or > 1) as we now show. 

2.3.2. For a > 1 the solution a.s. never reaches 

We now deal with the case a > 1 and show that in contrast to the case a < 1, even when 
5q is stable in probability or almost surely exponentially stable, the solutions of the SDE almost 
surely never reaches zero. 

Theorem 3. Let a > 1. Then for any xq g R we have 

F[x t * Oon t > 0] = 1 

That is, almost any sample path starting from a non-zero initial condition will never reach zero. 
Remark. A similar result (for a = I) in proved for instance in M8\ lemma 3.2]. 

Proof. Denote by t the first hitting time of zero of the process with initial condition x + at 
t — 0. Assuming that the conclusion is false, we can find a pair of constants T > and 9 > 1 
sufficiently large such that P(B) > where B is defined as: 

B = [t < T and \x(t)\ < 6 - 1 for all t < t) 

Moreover, under the assumption a > 1, both the drift and diffusion functions are locally Lipschitz- 
continuous, hence there exists a positive constant Kg such that: 

\Ax + g(x)\ V y(x) < K s \x\ for all \x\ < 6, t <T. 

Let now V(x) := |jt| _1 . We have for <\x\ <6 and t < T 

-CV < Kg(l + Kg)V(x) 

Let us now define r £ = inf{f > 0; |jc#| t (s, 0)}. This is a stopping time, so by Ito's formula and 
the optimal sampling theorem with the bounded stopping time r E A T, we have: 

■E[ e - K ><- 1 + K °X r > AT )v(x TtAT y\ = V(x )+ 

+ E[ J o e- K ^ +K " )s [-K e (l + Kg)V(x s ) + (£V)(x s )] ds] 
< \xo\~ 1 
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We know that for any oj e B, the time t e is smaller than T and moreover that x Tc — s by definition 
of B, hence we have: 



\xo\~ 1 > ne- Ke(l+Ke)(TcAT) V(^ T )] > M[e- K °V +K ° )T s- l l B ] 

P(B) < etar 1 **" 1 ^ 

This inequality is true whatever e (we recall that the definition of B did not involve s), hence 
valid in particular when e is arbitrarily small, and necessarily P(B) = which contradicts the 
definition of B. □ 

2.4. Stationary and quasi-stationary solutions 

We are now interested in identifying possible stationary solutions. We already mentioned 
that the null stochastic process (distributed as a Dirac delta measure at 0) is always a stationary 
solution of the system. Moreover, for a < 1, we have shown that the process is almost surely 
absorbed at zero in finite time. This property prevents from the existence of any stationary 
solution, and as far as the permanent dynamics is concerned, it is perfectly characterized by 
the absorption property. In order to further describe the behavior of the solutions, one may 
wonder how the solutions that did not yet reach zero are distributed. Stationary solutions of the 
process conditioned on never hitting zero are called quasi-stationary distributions. In the case 
■y(x) = o~\x\ a with a < 3/4, we can show that there exists a quasi-stationary distribution that 
attracts exponentially fast the conditioned process. In detail, we have: 

Proposition 4. Let us assume that v > 0, g is differentiate, g' is bounded at zero and smaller 
than Cx 2 P for some C € R at infinity, that assumptions (H1)\(H6) hold and a e [1/2, 3/4). The 



process conditioned on never reaching zero and with a bounded initial distribution converges 
exponentially fast towards a probability measure p ± for any positive (resp. negative) initial 
condition. These are called quasi-stationary solutions, or Yaglom limit. 

Proof. We apply a result due to Cattiaux and collaborators [ 15 ] in the domain of diffusion models 
in population dynamics. In their article, they prove consider a diffusion of type: 

dX t = dB,-q(X,)dt 

under a few conditions on q. Our problem can be brought to their framework changing variables 
and considering: X, = x)~ a /(cr(l - a)). Ito's formula allows putting our equations in that form, 
with 

q(x) = --Z- a f(Z) + —?— 1 - 
cr 2( 1 - a) x 

and where f(x) = Ax + g(x) and Z = (<x(l - a)x) ' ( - 1 ~ a \ Under the conditions that (i) the first hit- 
ting time of zero is almost surely finite and (ii) that the conditions: A = - inf 3 , e (0 >O o) q 2 (y) — q'{y) < 
oo and lim^oo q 2 (y) - q'(y) - °°, they show that the spectrum of the Kolmogorov backward 
equation (for the related operator £) has a purely discrete spectrum with no zero eigenvalue. 
Moreover, they show that a rescaled version of the related eigenfunction is integrable and expo- 
nentially attract all initial conditions under the conditions: 



f 

Jo 



- _GW dy < oo (6) 



q(y) 2 - q'(y) + A + 2 
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with G(y) = 2q(y) dy, 



( e _G6,) fify < oo , J xe- G(x)/2 dx < oo and ( e G(y) I e G(z) dzdy < oo (7) 

Jl Jo Jl Jv 

In order to apply their result, we demonstrate that all these conditions are satisfied in our 
framework. First of all, under the condition a < 1 and v > 0, we showed that the process x, 
solution of ([TJ almost surely reach zero in finite time (proposition^. This is hence also the case 
of X, since (1 - a) > in that case. Moreover, the derivative of the function q(x) reads: 

q'(x) = aZ- 1 f(Z)-f'(Z)- " ' 



2(1 - a) x 2 

Around zero, because of our assumptions on g, we have: 

2(1 - a) x 2(1 - a) x 1 

and hence 

2, x f f x ar(2-a) 1 

which tends to infinity as * — > + . At * = oo, we have: 

q \x)- q '(x)>^- 1 (l-s)Z^-^ 

for any s > and hence tends to infinity when x — > oo. In particular, A = - inf ve (o, M ) q 2 (y) - 
q'(y) < oo 

Let G(x) = 2q(y)dy. We now show that the condition (|6) is satisfied. At infinity, q(x) 
is greater or equal to c' x l+/3 ^ 1 ^ with C > 0. It is then easy to show that lim c ^ M = oo. 
Moreover, we have 

A := lim (G(x) log(x)) e (-oo, oo). 

Therefore, near zero, the integrand involved in the expression of Q(x) is proportional to y 2 ^. 
Under the condition or < |, condition (|6) holds. Showing that condition (|7} is valid is straight- 
forward at this point: G grows at least linearly at oo ensuring that e~ G ^dy < oo. Moreover, 
around zero, the function xe~ G ^ x ^ 2 is of order jt 1 ^^ which is integrable at zero for a < |, 
hence in particular for a < I. The main results of lfT31l (namely their theorem 5.2, proposition 
5.5, corollaries 6.1 and 6.2) apply, and ensure the existence and uniqueness of a quasi-stationary 
distribution attracting exponentially fast the law of the diffusion conditioned on not hitting zero 
with initial laws having bounded support. □ 

Let us now deal with the general case with diffusion coefficient y(x) and a > 1. In that 
case, the process is almost surely never absorbed at zero. Possible stationary solutions with 
smooth probability density function p with respect to Lebesgue's measure satisfy the forward 
Kolmogorov (or Fokker-Planck) equation: 



^ ((Ax + g(x))p(x) -~ (y(x) 2 p(x)j^ = 0. (8) 
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i.e. 

1 2 , 

(Ax + g(x) - y (x)y(x))p(x) - -y(x) p (x) = K 

for some K e R. If we can find solutions to these equations subject to the conditions p(x) > for 
all x e R and L p(x) dx — 1 (acceptable solutions), then these are probability density functions 
of stationary solutions of the diffusion equation. For simplicity, we consider that assumption |(H5)| 
is valid, i.e. y(x) only vanishes at x = 0. In that case, since do is a solution to the SDE ([T}, 
stationary solutions will either have support on R + or R~. In the case where the dynamics 
has additional singular points, the Kolmogorov equation will define solutions in each interval 
delimited by the singular points (see section [33) . The Kolmogorov equation is a linear ordinary 
differential equation with constant source term. The homogeneous solutions (with K = 0) is 
given by: 

2-tv , . g(y) 



/,„(,-) = Cy(.vr : «p( J (^+2^),/v 



defined for x + 0. Particular solutions are given by: 
pM - (-2K J "exp(-2 1' |1 + 2^ *)* + r(l) 2 p(l) 



■or 



We have the following: 
Proposition 5. In addition to the stationary distribution 5q, we have: 
• For a — I and: 

- A < a 2 /2 there exists no integrable solution of the Kolmogorov equation 

— A > cr 2 /2, under assumption (H6) there exist two integrable solutions of the Kol- 



mogorov equation given by Z + po(x)l x> o and Z Po(x)l x <0 where Z + and Z generi- 
cally denote normalization constants. 

9 6 I I 

— A > cr /2, v < and 6 < 1 + j (case where assumption (H6) is not satisfied) there 



does not exist any stationary solution (we saw that the solutions a.s. blow up infinite 
time) 

For a > 1 and: 

— A < the Kolmogorov equation has no integrable solution 

— A > and under assumption \(H6)\ the Kolmogorov equation has two integrable 
solution given by Z + pq(x)1 x> q and Z pq(x)1 x< q. 

— A > 0, v < and 5 < 1 + f (i.e. when (H6) is not satisfied), the Kolmogorov equation 



has no integrable solution. 

REMARK. In the case a < 1, we already mentioned that no stationary solution can exist. We can show 
using the expressions of the solutions to the Kolmogorov equation that no solution is integrable. 
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Proof. For a > 1, let us analyze the integrability of po(x). Around 0, the function x is not 
integrable. The leading term of the exponential around zero is given by ■^— ) x 2{x ' a \ hence 
diverges at zero. If A < 0, this term tends to +oo and is not integrable, and if A > 0, the 
exponential term tends to zero and the function is hence integrable at and p(Q) = 0. For 
x — > oo, we distinguish two cases: 

• for 5 < 1 + /J/2 and v > 0, the leading term of the exponential in the expression of po(x) 
is given by -vz 1+/3 ~ 2<s /(l + /3 - 26) which tends to -oo for v > fast enough ensuring 
integrability of po(x) at oo. For v < 0, this term tends to oo and po(x) is not integrable at 
infinity. 

• when v < 0, the leading term of the exponential at infinity is of order —vz x+ ^ 2s /(l +f3—28) 
which tends to +oo and hence the solution to the Kolmogorov equation is not integrable at 
infinity. 

• For 5 > 1 + (3 /2, the terms inside the exponential tend to zero, and the integrability relies 
on the integrability property of x~ 2s , always integrable at infinity. 

We hence proved that acceptable solutions related to K = exist for A > and no solution 
existed when A < 0. Moreover, there is no possible choice of K + that can overcome the 
exponential divergence at zero. Indeed, the only possible choice of K would be the only constant 
preventing divergence at zero, i.e.: 

However, the integral term of the righthand side diverges, hence the above relationship cannot be 
satisfied and there is no acceptable solution for A < 0. Eventually let us remark that for A > 0, 
any choice of K prevents integrability at infinity. 

For a - 1, we again analyze the integrability properties of po(x). Around x — 0, p$ behaves 

like exp(2 J* dyj. The term in the exponential behaves like ^ky". The integrability at 
zero hence depends on the exponent -2(1 - -^): po is integrable if the exponent is strictly larger 
than -1, i.e. 2A/cr 2 > 1. In the case 2Ajcr 2 < 1, the non-integrability is due to a polynomial 
divergence, and hence might be compensated by a suitable choice of K. The only possible choice 
of K compensating the divergence at zero would correspond to: 

and here again the integral diverges at zero, and hence there is no such solution. 

At infinity, the exact same analysis as done in the case a > 1 applies, and hence we obtain 
that po is integrable at oo under assumption |(H6)| □ 

2.5. Heuristic discussion 

We hence proved that as a function of the value of the coefficient a, the system can be in 
one of two substantially different regimes: in the case a < 1 the singular point is always stable, 
and reached in finite time, and for a > 1 the stability of the singular point is not affected by the 
presence of noise and only depends on the stability of the fixed point in the noiseless system. 
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The case a — 1 constitutes the transition between these two regimes, and in that case the stability 
of the singular point depends both on the eigenvalue of the Jacobian matrix of the drift at the 
singular point and on the level of noise. 

In order to understand heuristically these stabilization and destabilization phenomena, let us 
focus on the behavior of the system in the neighborhood of the singular point. In the case a > 1, 
the diffusion coefficient vanishes faster than the drift coefficient at the singular point, and the 
system locally behaves as if there were no noise in the system. In particular, the stability of the 
singular point remains unchanged from that of the noiseless system. On the contrary, when a < 1, 
the diffusion coefficient vanishes slower than the drift. When the system approaches the singular 
point, it is hence mainly driven by random fluctuations related to the diffusion coefficient. When 
the system is brought away from the singular point, noise increases, and when it is brought 
towards the singular point noise decreases. This induces a ratchet-like phenomenon: consider to 
points A and B close of the singular point, A being closer than B. The time it takes to go from A 
to B with power diffusion coefficient is way larger than the time needed to go from B to A. This 
behavior of the noise close of zero explains the stability in probability, and acts as a weak drift by 
'stabilizing' the solutions when they approach the singular point. The intermediate case a = 1 is 
precisely the transition where these two effects compensate. When the noise coefficient is large 
enough compared to the eigenvalue of the Jacobian matrix at the singular point, the noise effects 
described for a < 1 dominate and stabilize the fixed point, and on the contrary the deterministic 
phenomena govern the dynamics similarly to the case a > 1. 



3. The stochastic pitchfork, Hopf and saddle-node bifurcations 

Now that we characterized in detail the solutions of equation ([TJ, we are in a position address 
the dynamics of three prominent bifurcations occurring in the study of dynamical systems: the 
pitchfork, Hopf and saddle-node (fold) bifurcation with power diffusion coefficients. 

3.1. Analysis of the stochastic pitchfork bifurcation 

In this section we first study the dynamics of an SDE with drift given by the normal form of 
the pitchfork bifurcation, before addressing the universality of these behaviors. 

3.1.1. Dynamics 

We consider the supercritical stochastic pitchfork bifurcation with power diffusion. 

dx, = (Ax, - jc, 3 ) dt + o-\x,\ a dW t (10) 

for x, e R, with initial condition Xq at t = 0. In this equation, A is a real parameter, cr is a 
non-negative parameter and a > I. In the case cr = 0, it is well known that the solution x = 
is a stable fixed point for any A < 0, and unstable for A > 0, and two additional stable equilibria 
+ V/l exist in the region A > (see e.g. 0]). 

For cr + 0, the model clearly satisfies assumptions (H1)||(H6) with v = = 1 > 0, k - 



B — 2, d — cr and 6 = a. Moreover, it satisfies the additional assumptions used in proposition [4] 
since g(x) = -x 3 is differentiable with bounded derivative at zero and divergence at infinity 
upperbounded by x A . All results of section [5] apply and we have: 

Theorem 4. • For a < 1, the fixed point Sq is stable in probability. Moreover, any solution 
almost surely reaches zero in finite time. Quasi-stationary distributions exists for a e 
(1/2,3/4). 
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• For a — 1 and : 

— A < cr 2 12, do is almost surely exponentially stable and no additional stationary solu- 
tion exist. 

— A > a 2 /2, 6q is unstable in probability, and two stationary solutions stable in the first 
approximation exist. 

• for a > 1 and: 

— A < 0, <5o is stable in probability and no other stationary solution exist. 

— A > 0, do unstable in probability and stationary solutions exist. 

Proof. Most results presented in the theorem are direct applications of the results demonstrated 
in section [2] The only assertion to be demonstrated is the stability in the first approximation of 
the additional stationary solutions. 

In the case a = 1, these solutions have the expression: 

where we made explicit the normalization constant. The linearized equation reads: 

dv, = (A-3 X 2 )v, dt + o-v, dW, (11) 



which is integrated as 



cr 2 f 

v, = exp [(A - —)t - 3X 2 ds + o-W,]v , 
2 Jo 

and hence the Lyapunov exponent is given by: 

limsup - log |v,| = (A - — ) + o--^ -- I 3 X 2 S ds + - log(v ). 

>oo t L t t Jo t 

Using the strong law of large numbers ifTTl Chap.IV.5], the term - t j Q 3X 2 ds converges a.s. to- 
wards 3E[X 2 ] which can be computed in closed form, yielding the following expression of the 
Lyapunov exponent: 



1 ' io- 1 ' 1 



2 2 

cr cr 
= (A- T )-3(A- T ) 

cr 2 

= -2(A-—)<0 

which ensures stability in the first approximation of the new stationary solutions since A > 

□ 
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(a) or = 0.6 (b) a = 0.85 



Figure 2: Distribution (ordinate truncated for legibility )of the trajectories that did not reach zero in the subcritical pitch- 
fork model with A = -0.5, o~ = 0.5 different values of a, computed at 6 times. (a)a = 0.6: stabilization on a quasi- 
stationary distribution that do not depend on time, a = 0.85: the repartition varies as a function of time and all trajectories 
approach zero, arguing for the non-existence of quasi-stationary distribution. Simulations were done using a particle 
method with 100 000 particles, using the Euler-Maruyama scheme with dt = 0.01. 



Let us now further describe the dynamics of the solutions as a function of a. 

• Case a = 1: We observe that for any A > 0, the unstable deterministic fixed point becomes 
asymptotically exponentially stable when the noise parameter cr is large enough. For A > 

two symmetrical stationary distributions appear. For ^- < A < a 2 the stationary 
distribution concentrates at zero and has a non increasing density diverging at zero. For 
A > cr 2 , the probability density function of the stationary distribution vanishes at zero 
and has a unique maximum reached for x — + V/i - cr 2 (see Figure lij. There is hence a 
qualitative transition at A = cr 2 , or P-bifurcation. In comparison with the deterministic 
bifurcation, the loss of stability is delayed and noise tends to stabilize the saddle point. 

• Case a > 1: the deterministic picture is qualitatively and quantitatively recovered: for 
A < 0, 6o is stable in probability and is the unique stationary solution, and for A > 0, 6q 
is unstable in probability, two additional stationary solutions appear, presumably stable as 
shown for a — 2. The stationary probability density functions for A > vanish at zero and 
reach a maximum at X„, the solution of 



(<r 2 a-AX 2 - 2a +X- 2a+4 ) = 0. 
The behavior of x m when A is close to zero follows one of the two regimes: 

/ . > 1 1(2(1-2) , , , 

- for a < 2, X m ~ 1^ J . For instance for a = | we get X m = -\cr 
\ V9cr 4 + \bA 

- for a > 2, X m ~ V/i and is insensitive to noise parameters cr and a. For a — 3 we 

1 V-1+ VT+i2ir 2 ^ 



obtain the explicit form X m - 



- for a — 2, this maximum can be computed explicitly, and is reached for x = J 



l+2o- 2 



Surprisingly, in the limit A — > there is a discontinuity in the behavior of X m as a function 
of a at a — 2. Nevertheless, when A > this discontinuity is smoothed out. We conclude 
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Figure 3: (a) Codimension 2 stochastic bifurcation diagram of the pitchfork bifurcation: orange: zero is stable, blue: 
zero is unstable, yellow: zero is stable and is reached in finite time, (b) Sample path trajectories in all 6 regimes for the 
stochastic pitchfork bifurcation. Parameters (A, a, cr): (bl) (10, 0.6, 1) (quasi-stationary distribution plotted in Fig.|2|a)) 
(b2): (2, 1,1.2), (b3): (1, 1.5, 1.2), (b4): (0.5,1,1.2), (b5): (-1,1.5.1.2). In (b2) and (b3) the peak of the stationary 
distribution is depicted in dashed black, (c) Representation of stationary densities for the stochastic pitchfork system 
with linear noise, as a function of the parameter A. Red square: 6q is stable. Blue and colored surfaces display the 
probability distribution p,\(x) of the stable stationary distributions, and the black line represents the maximal value of the 
distribution. The case a = 1 was also studied using the theory of random dynamical systems in (9)- 



that in the case a > 2, the behavior of the system is close from the one of the deterministic 
pitchfork bifurcation and the stationary distribution escapes zero as a square root of the 
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parameter. But when 1 < a < 2, this is no more the case: the scaling of distribution 
peak behaves as A}^ 2 " -2 ^, hence mucn slower than VX As cr is increased, the peak of the 
stationary distribution gets closer to zero. 

• Case a < 1: the solution 5q is always stable in probability. This result can appear relatively 
surprising at first sight. Indeed, in the deterministic case, A is the exponential rate of diver- 
gence from the solution 0. However, adding a (possibly small) diffusion term proportional 
to \x\ a with a < 1 stabilizes 6q in probability whatever the value of the noise intensity 
cr + 0. This observation, added to the fact that the solution is not exponentially asymp- 
totically stable (though stable in probability) raises the question of how this convergence 
occurs. Figure|3jbl) presents a typical sample path of the process. Starting from a positive 
initial condition, we observe that the solution is evolving in the half-plane x > and does 
not show any extinction clue. However, it suddenly reaches zero where it is absorbed after 
that random transient phase. This perfectly illustrates a typical behavior of the solutions of 
the pitchfork equation for a < 1: the extinction time was shown to be almost surely finite. 
Moreover, trajectories that did not hit zero are distributed according to a quasi-stationary 
distributions as long as a < 3/4. 

Remark. Let us briefly discuss the case of the subcritical stochastic pitchfork bifurcation: 

dx, = (Ax, + x]) dt + cr\x,\ a dW,. 

The analysis of this system can be treated in the same fashion as the supercritical case. In the deterministic 
case (cr = 0), is stable for A < and unstable for A > 0. Two additional unstable fixed points ± y—A 
exist in the region A < 0, and for A > the solution of the deterministic pitchfork equations blow up in 
finite time. In the case a = 1, the analysis of the stochastic case gives the following conclusions. The fixed 

2 

point of the subcritical stochastic pitchfork equation is almost surely exponentially unstable for A > S T , 

2 J 

asymptotically unstable in probability for A > — ^- and asymptotically stable in probability if A < — No 
extra stationary solution appear in that case. 

3.1.2. Universality properties 

Similarly to the analysis of nonlinear dynamical systems, the analysis of the pitchfork equa- 
tion reveals universal features of a wide class of systems. We formally show the universality 
property of the pitchfork bifurcation for general SDE on R: 

dx t — f(x t , A) dt + g(x t )dW, 

with a smooth / (at least three times continuously differentiable). We consider the dynamics of 
the system for parameter values A in a neighborhood of Aq e R. Let us assume that: 

• there exists v > and k > such that for all A in a neighborhood of Aq, f(x, A) < -vx l+K , 

• for any A, f is an odd function: f(-x, A) = -f(x, A) 

• g((Uo)*0 

Defining £ = (3(A)x, with /3(A) = J i 1^(0, A)|, we have: 

dtt = (^(0, X£ t + + V(g„ A) J dt + f3(A)g dW t . 
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where ^(x, A) = 0(x 5 ) and s = sign(^-(0, A)) which is well defined around A = Ao. 

Let us now assume that g{x) ~ cr\x\ a around x = and g(x) > for x + 0. We can hence 
apply the results developed in the previous section to characterize the dynamics of this general 
system. Pitchfork-like transitions occur at A — A§ when: 

• a > 1, %(0,Ao) = and gj(0, Ao) * or 

• a = 1, f ((Uo) = f and ^(0,^,) * 0. 



We now put in good use this reduction in the following example. 



Example 1. Let us consider the Wilson and Cowan firing- rate model [19], classically describing 
the behavior of a cortical column. In that model, the voltage v of a typical neuron in tlie column 
is solution of the equation: 

dv, = iy--v, + JS(v,)j dt + y(v,)dW, 

where S(x) is a smooth sigmoidal function describing the voltage-to-rate transformation, and J 
is the typical connectivity strength. We choose here S (x) — tanh(g x) with g > the sharpness 
coefficient of the sigmoid and y(x) - o~\x\ a with a > j. This equation can be locally reduced to 
the normal form of the pitchfork bifurcation since the flow is symmetrical, and we obtain using 
the same transformation as performed above: 

dt, = {i-\ + ~ £ + m /,*)) dt + 2 1 -V|f t |W, 

For a < 1, the fixed point is almost surely reached infinite time, for a > 1, the fixed point is 
stable in probability for gjr < 1 and unstable in probability otherwise. For a — I, the fixed point 
is stable in probability if — 1/r + gj < cr 2 /2 and unstable otherwise. This behavior is illustrated 
numerically in Figure^ 

3.2. The stochastic Hopf bifurcation 

The Hopf bifurcation for smooth dynamical system is topologically equivalent to the two 
dimensional normal form: 

dx 77 dy 77 

— - Bx-y + sx(x +y ) ; — - By + x + ey(x + y ) 
At At 

where /3 is the bifurcation parameter and s = ±1 governs the type of bifurcation. It can be 
conveniently written in complex form for z, = x, + iy, with i 2 = -1: 

% =08 + i + e |z| 2 )z. 
At 

The null solution (x — Q,y = 0) is always solution of the equations. It is exponentially stable 
for any /3 < and exponentially unstable for any /? > 0. If s = — 1, the system presents stable 
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(a) a = 0.8 




(b) a = 1 



(c) a = 1.5 



Figure 4: Wilson and Cowan system with power diffusion function: empirical mean m of Vf at time T = 200 for 1 000 
independent sample path as a function of g. a < 1: the fixed point is reached in finite time whatever g, a = 1: zero 
looses stability for g = 1.5 and a > 1: zero looses stability for g = 1. Due to the finite simulation time, the transitions 
appear smooth in these simulations. Euler Maruyama scheme with time step dt = 0.001, 100 values for parameter g. 



oscillations for/? > and the bifurcation is termed supercritical, and if s = 1, unstable periodic 
orbit exist for (3 < and the bifurcation is termed subcritical. 

We study in this section the stochastic Hopf bifurcation with multiplicative noise. For sim- 
plicity, we will consider that the stochastic perturbations are driven by a single Brownian motion 
W,. Let <x* = cr + ifi be a noise parameter, the stochastic Hopf normal form equation reads (in 
complex form): 

dZ t = (/3+i + s\Z,\ 2 )Z, dt + cr*\Z,\ a - l Z, dW„ 



which corresponds for Z, - X, + \Y t to the bidimensional system: 

dX, =(J3X t -Y t + sX, (Xf + if)) dt + (o-X,-fi Y,)\Z,\"- 1 dW, 
dY, =(/3Y,+X,+ sY, (Xf + Y?)) dt + (^X, + cr Y,)\Z,\ a - y dW, 



(12) 



where and s correspond to the parameters of the Hopf bifurcation and fi and cr are two param- 
eters globally governing the amplitude of the stochastic perturbation. 

The deterministic process Z, = corresponding to X, = 0, Y, — for all f > is solution 
of ( p"2"} and is univocally defined by the fact that the modulus of Z, is null. We denote by R t the 
modulus of Z, and by 0, its argumenj^] 

Lemma 1. The modulus of the variable R, — \Z,\ and the argument 6 t satisfy the equations: 

idR, = {BR, + f Rf*-* + sR^ dt + o-Rf dW, 
\d9, = (-o-iuRf '- 2 + l)dt+ nRf- 1 dW, 

Proof. Let (X t , Y,) be a solution of the Hopf equations ( fT2| ). We apply Ito formula to the variable 



2 R, = \Z t \ = ijxf + Yj , 8, given by the arctan(^) when X t > 0, added or subtracted n when X, < depending on 
the sign of Y,, and ±n/2 if X, = 0, ±¥, > 
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R,= Jx] + Y 2 : 

dR, = —(X,dX, + Y,dY,) + — (Y 2 d(X), + X 2 d(Y), + 2X, Y,d{X, Y) t ) 
Rt 2R* v ' 

= {^{pX 2 -X,Y, + eX 2 (X 2 + Y?) + /3Y? + X, Y, + eY 2 (X 1 , + if)) 

+ -j-^\Yj{crX t - p Y,) 2 + X 2 (a- Y t +n X,) 2 - 2X, Y, (aX t -fx Y,) (crY.+p X,)]} 



dW r 



+ -L- \X, (aX t -fx Y,) + Y t (crY t +fx X,)} 
i(0R? + sRt) + ^ 2 RT 2a - 2 ) dt + <W 

= + + S^j (if + cr/?? t/W, 

The argument 9, is given by 6, = arctan(F r /X,) plus possibly constants depending on the sign of 
X, and Y,. Applying Ito's formula again yields: 

dB, = (1 - o-pR 2{a ~ l) )dt + pR^ l dW, 

which ends the proof of the lemma. □ 

It is important to note that the equation on the modulus is uncoupled of the phase equation 
on 0,. The modulus of (X,, Y,) is therefore solution of a stochastic differential equation of type 
bifurcation (fill. Moreover, when a = 1, the equations take the simpler form: 

{dR, = {(J3 + lj)R, + sR^ dt + crR t dW, 
\d6, = (l-o-p)dt+pdW, 

and hence the variable R, is solution of a stochastic pitchfork bifurcation. By application of 
theorems |4] and remark [3. 1.1 1 it is easy to establish the following: 

Theorem 5. For a — I, the null solution of the supercritical Hopf equations is almost surely 

2 2 ^2 2 

exponentially stable if fj < Z and asymptotically stochastically unstable if fj > ^ . ^ n tnat 
case, there exists a new stochastically stable stationary solution with distribution: 

The null solution of the subcritical Hopf equations is almost surely exponentially unstable if 
P > 2 M . It is asymptotically stochastically unstable if [3 > — V 1 and stochastically stable if 

For <x + 0, the modulus R, is a.s. exponentially attracted towards when /? < ~^ , There- 
fore, as the imaginary part of the noise, cr, is increased, the fixed point gains stability. For 
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Figure 5: Simulations of sample path trajectories of (X,, Y t ) for the stochastic Hopf bifurcation with a = 1, jx = and 

2 

<r = 1 for different values of ft (deterministic phase case). (A):/? = 0.3 < zero is a.s. exponentially stable; (Al): 
trajectories (blue: X t , red: Y t ), (A2): phase plane, (A3): phase 9, is deterministic equal to V. the solutions rotate with 
pulsation one. This is true of all three cases (A), (B) and (C). (B): fi = 0.7 e [^p, <r 2 ]: Zero is unstable and the stationary 
modulus distribution is peaked at zero; (Bl): trajectories, (B2): phase plane, (B3): distribution of the modulus R t . (C): 
P = 4 > o" 2 : zero is unstable and the distribution of the modulus has a peak around R* > (C3): the trajectories display 
perfectly periodic oscillations (C 1 , C2) with random modulus, a highly regular behavior for a stochastic system. 

P e (^-, ^-^-), the fixed point zero is unstable and solutions do not converge towards the 
fixed point zero. The distribution of the moduli R, converge towards the stationary distribu- 
tion Pr(x) heavily charging zero, and solutions stay close from the fixed point with random 
stochastic variations. As soon as /3 > — ^ , the solutions have their modulus centered at 

R* := - ^"2^- and present stochastic oscillations with a deterministic pulsation (-1 + cr/i) 
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and random phase given by a Brownian motion scaled by y. When y — 0, the phase is hence de- 

2^2 2 ^j-2 2 

terministic and the system rotates with pulsation 1. When /3 e ( J 1 , - ), the stationary mod- 
ulus of the solution has a distribution peaked at zero and no actual oscillation can be observed. 

When p > 2<T ^ > tne stationary modulus of the solution has a maximum at R* - - ^-f^- 
and the distribution is null at zero: solutions hence correspond to perfectly periodic oscilla- 
tions with random amplitude. Interestingly, when cry crosses 1, the direction of the oscillation 
changes: for cry < 1, the oscillation is counter clockwise, and for cry > 0, the oscillation is 
clockwise. This is a new type of bifurcation due to noise that we term orientation bifurcation. 
This phenomenon is illustrated in Figure|6] In the case cry = 1, the pulsation is null and the phase 
is given by yW,: it has a probability 1 /2 to rotate clockwise and 1/2 to rotate counterclockwise, 
and changes the sense of rotation. The rotation number is equal to zero. 




(a) cr = 0.5 (b) cr= 1 (c) 0-= 1.5 



Figure 6: Noise-induced inversion of the rotation pulsation: a = i, ft = 4, n = 1 for different values of cr. Blue: phase 
(6 t ) as a function of time and black: deterministic pulsation 1 - cr/i. (a) counter-clockwise rotation, (b): no rotation, 
trajectories appear stochastic and the phase is a Brownian motion, and (c): clockwise rotation. 

We observe in the purely imaginary noise case (cr = 0), R t satisfies the deterministic pitchfork 

bifurcation equation. For /3 > y , R, - + y is the only attractive solution: the solution of 
the stochastic Hopf bifurcation asymptotically lives on a deterministic circle. Its rotates on this 
circle with a stochastic phase given by a Brownian motion multiplied by y and a deterministic 
pulsation equal to 1 (see Figure|7]i. 

Let us eventually note the opposite effects of the real and imaginary values of the noise 
parameter: for fixed ft, increasing cr tends to stabilize 0, whereas increasing y has the effect of 
destabilizing 0. These two effects compensate each other, and when choosing cr — y, the loss 
of stability of zero arises at ft — precisely as it is the case of the deterministic bifurcation. 
However, cycles appear for f3 > cr 2 /2 in that case, and hence we observe a systematic delay in 
the apparition of cycles for cr + 0. 

Beyond these interesting effects of noise, we want to underline the fact that for a — 1 the 
global picture of the deterministic Hopf bifurcation is nevertheless conserved, since as /3 is in- 
creased, a limit cycle appears, around which the solutions stochastically rotate. The picture will 
be relatively different in the nonlinear noise case a + 1 . In that case we prove the following: 

Proposition 6. In the case of the supercritical (e — —I) stochastic Hopf bifurcation, we have: 

• For a < 1, the fixed point is: 

- stable in probability for \y\ < \cr\ 

- unstable in probability for \y\ > \cr\ 
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Figure 7: Simulations of sample path trajectories of (X,, Y,) for the stochastic Hopf bifurcation with a = 1, cr = and 

fi = 2 (deterministic modulus case) for different values of fl. (A):fi = -2 < — ^j-: zero is a.s. exponentially stable; (Al): 
trajectories (blue: X t , red: Y t , black: R t s), (A2): phase plane, (A3): argument 8, in blue, deterministic drift in black. 
(B): = > — K- : the modulus converges to a constant value -y//i 2 /2 = V2; (Bl): trajectories, (B2): phase plane, (B3): 
Argument 0, (blue), deterministic drift (black). 



Moreover, almost all sample path reaches zero infinite time. 
• For a > 1, the fixed point is: 

— stable in probability for /3 < -y 

- unstable in probability otherwise. Two additional stationary solutions exist in that 
case. 

The solutions almost surely never reach zero. 

Remark. Let us emphasize the fact that the case a < 1 is substantially different from the general case: 
indeed, we had seen that the singular point was always stochastically stable. Surprisingly here the stability 
of zero depends on the relative values of the real and imaginary parts of the noise. 

Proof. The case a > 1 falls in the general analysis developed in section|2] Indeed, the d iffusio n is 

and 



of the form of equation ([TJ with g(x) = yx 2ff 1 - x 3 which clearly satisfies assumptions 



(HI) 



(H2) The general analysis hence applies and directly leads to the conclusion of the proposition. 
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In the case a < 1, the equation was not covered in the general analysis, since the drift includes 
a term yx 2 "' -1 and 2a - 1 < 1. However, the same principles may be applied to that case. First 
of all, existence and uniqueness of solutions is proved using the same arguments, and the a.s. 
boundedness of the process as well. Hence there exist a unique solution to the stochastic Hopf 
bifurcation process with a < 1, which is almost surely bounded for all times. The stability of 
zero is investigated as follows. For fj. > cr, the Lyapunov function V(x) = - log(x) is C ( 2 around 
zero, positive, diverges at zero, and 

x,/w = K + y) A2+0(A2) 



and is hence negative around zero, proving that zero is unstable in probability for /i > cr. 

2a- 1 



For /i < cr, we define the Lyapunov function V(x) = x* with <1> = -= § > 0. This function is 



Cy and we have: 

£V(x) = $ fcj + (* - l)y ) x 2a - 1+m + o(x 2a - 2 ^) 

and hence is negative close of zero. Using the same techniques as in the general theory, it is easy 
to conclude that the singular point is unstable in probability. In order to prove that the solutions 
almost surely hit zero in finite time, we proceed as in the proof of proposition^ We know that the 
process is almost surely bounded, and analyze Feller's scale functions. It is clear that p(x) — > oo 
as x — > oo following the proof of proposition]^ the proof only uses the higher order term -x 3 . 
We hence need to ensure that lim x _ > o+ v(x) < oo. This property is related to fine properties of the 
scale functions close from zero, which is now governed by the term /i 2 /2 x 2 " -1 . We hence have, 
around 0: 

2 ^ x 2 ~ 2ff x 2(2-a) 

2<t z cr 2(1 - or) 2(2 - a)cr z 

and therefore we can obtain the explicit expression of p(x): 



„4-2o- 



More importantly, in order to characterize v, we consider the term G(y)-G(z) which is equivalent 

JL 

2a- 2 

,2 



close from to: G(y) - G(z) ~o+ yi log(-). This implies that: 



J exp(2G(y)-G(z))dz~o+ J exp^log(^)J* 

1 //y\^ 2/cr ' 



1 -n 2 /cr 2 



- 1 



and hence: 



v(jc) 



cr 2 \fl 2 /cr 2 + 1 

This function will hence remain bounded as x — > oo, ensuring using Feller's test for explosion that 
the first exit time of the interval (0, oo) is almost surely bounded. The almost sure boundedness 
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of the process allows concluding that the probability first hitting time of zero is almost surely 
finite, which ends the proof. □ 



Therefore, similarly to the case of the stochastic pitchfork bifurcation, we observed that for 
a > 1, the noise does not affect the stability of and the loss of stability precisely appears for 
/3 = as in the deterministic case. The behavior of the system beyond the bifurcation point 
depends on the shape of the stationary distribution, and in that case again the solution reaches a 
maximum for a positive value, reproducing qualitatively the behavior of the deterministic pitch- 
fork bifurcation. 

For a < 1, the almost sure absorption at zero of the solution persists. However, the stochastic 
stability of the fixed point zero now depends on the relative values of p and <x. This distinction 
is not fundamental as long as we are concerned with sample path properties, and similarly to the 
general case, almost any sample path will eventually be absorbed at zero in finite time. 

Let us eventually remark that the universality properties of this normal form are harder to ob- 
tain than in the case of the pitchfork case already treated and the saddle-node bifurcation to come. 
This is not only due to the fact that we are in higher dimensions, but also (and this difficulty is 
probably harder to overcome) because the generic reduction to normal form in the Hopf bifurca- 
tion involves a nonlinear change of variable, which is not the case of the pitchfork or saddle-node 
bifurcation (see [2|). Since we consider diffusion processes, the process obtained through non- 
linear change of variable will satisfy, using Ito formula, a SDE where the drift function involves 
a term related to the diffusion, and it is hence much more complex to obtain a suitable change of 
variable. Another approach would consist in showing that some function of the process satisfies 
an SDE of the type studied in section|2] Let us formally give the kind of computation that would 
be involved, considering a two-dimensional system: 

dX, = F(X,) dt + G(X,)dW, 

where G(x) is precisely of the complex noise form considered above. The modulus of the solution 
R, satisfies the SDE: 

dR, = i- (x]F\X,) + X?F 2 (X t ) + ^-fl,J dt + (XjG\X t ) + X 2 G 2 iX t ))dW, 

where X @\ F 1(2) and G 1(2) represent the first (second) component of the vector field. The ques- 
tion that now arises is whether this system can be reduced to a pitchfork equation with possibly 
higher order terms, which can be checked in a case-by-case basis, as we now do for our choice 
example. 

Example 2. Let us consider a Wilson and Cowan neural network composed of an excitatory and 
an inhibitory population. The network equations in that case read: 

dv) = (-v) + S{gv]) - S(gv 2 )) + {o-v) - pvf)\v,r l dW, 
dvj = (-v 2 + S(gvj) + S(gvj)) + (o-v 2 + pv})\v,\ a - l dW, 

with S(x) = tanh(gx) and \v t \ = -yj(v]) 2 + (Vj) 2 . In the deterministic case, it is straightforward 
to show that the system undergoes a Hopf bifurcation at g — 1. Let us now consider the equation 
satisfied by the modulus p, — \v t \. Using ltd' s formula we obtain: 

dp, = ((-1 + g) Pt + p 2 /2p 2a - 1 + g{v),V 2 )) dt + Orp^dW, 
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with 



g(x,y) = -fV + / + x 3 y + xy 3 ) + 0(p 4 ) 
3p 

= -^((x+y) 4 - 2x 2 y 2 ) + 0(p 4 ) <~P 3 + 0(p 4 ) 
3p 12 

The argument satisfies the equation: 

dO, = (g- o-w 2a - 2 + 0(p 2 )) dt + pp'^dW, 
The dynamics of the system can hence be analyzed using the above analysis: 

• a = 1, zero is stable in probability when (-1 + g + p 2 12) < <x 2 /2 and unstable otherwise. 
At this point, a stochastic Hopf bifurcation arises, and oscillations appear, with pulsation 
locally equivalent to g — pa (with order two corrections in p) and stochastic phase given 
by a Brownian motion multiplied by p. The oscillation orientation bifurcation takes place 
forg =po-. 

• For a > 1 zero is stable if and only if(—l + g + p 2 72) < 0. 

• For a < 1, the fixed point is stable in probability if\p\ < \o~\ and unstable otherwise. The 
first hitting time of the singular point is almost surely finite in both cases. 

The rotation orientation transition is illustrated in Figure^ 

3.3. The stochastic saddle-node bifurcation 

The third prominent universal codimension one bifurcation of equilibria in deterministic 
ODEs is the saddle-node (or fold) bifurcation. It corresponds to the normal form i = -x 2 + a 
(sometimes considered x 2 - a which is equivalent to the normal form considered here through a 
change of time). This dynamical system has no equilibrium when a < (solutions blows up in 
finite time) and two equilibria for a > 0: -\/a which is stable and - -\/a which is unstable. The 
system has multiple singular points, and these depends on the bifurcation parameter. Choices of 
diffusion functions vanishing only at -\/a or at - -\/a enter the application domain of the general 
theory developed in section [2] Let us consider now a case, contrasting with the previous cases 
treated, where the diffusion function vanishes at both singular points: 

dX, = (-X 2 + a)dt + o-\X 2 - a\ a dW, (14) 

For a < 0, there are two stationary distributions corresponding to the equilibria -\/a and - -\fa. 
Denoting y, + = x, - yfa (vanishing at x = -\fa) and y~ — x, + yfa, it is easy to show that these 
variable satisfy the equations: 

(dyt = (-2 V5tf - (yt) 2 ) dt + cr\-2 ^ - (y?) 2 \ a dW, 
\dy; = (2 V5yr - (yt) 2 ) dt + cr\2 ^yt - (yt) 2 f dW, 



In both case, the drift and diffusion coefficients satisfy the assumptions (HI) (H2) (H3) 



(H4) 



and (H6) and obviously assumption (H5) is not valid. Assumptions (HI) (H2) and (H6) are 
trivial since the function corresponding to g(x) in the general cases is equal to -x 2 , hence p = 
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Figure 8: Dynamics of the Wilson and Cowan two-populations system {13\ . (a)-(d): vj (blue) and vj (red) as a function 
of time; (e) and (f): 8 t (green) and deterministic drift (black) as a function of time, (a) and (d) correspond to a < 1: for 
< \<t\, the singular point is stable (d) and unstable otherwise (a). The other figures correspond to a = 1 and illustrate 
the change of rotation orientation bifurcation at crfi = g. All the figures correspond to g = 1 . 



v = — 1 < and /3 — k — 1. The diffusion coefficients denoted y(x) in the general cas e satisf y 

is 



for both the equation on y + and on y : y(x) ~ x -^o* cr(2 ~\ /a) a x a , ensuring that assumption (H3) 
satisfied, and y(x) ~ A -^oo cr\y\ 2a ensuring that assumption (H4) is satisfied with d — cr and 5 -2a 



In the deterministic case, for a > and an initial condition smaller than - -\/a or for any 
initial condition when a < 0, the solutions blow up in finite time to -co, and in usual applications 
additional, higher order confining terms prevent the blow up. Interestingly, in the present case, 
the presence of noise can prevent blow up as we now show: 

Proposition 7. The solutions of the stochastic saddle-node equation ( |14| > are defined for all times 
or blow up infinite time under the conditions: 

• For j < a < 1 we have: 

— if a < the solutions almost surely blow up infinite time to — oo 

- if a > 0, the solutions with initial condition smaller than — yfa almost surely exit the 
interval (—oo, — yfa) infinite time. 

• for a > 1, the solutions almost surely never blow up, and for a > 0, the exit time of 
(—oo, — yfa) is almost surely infinite. 

Proof. We demonstrate this proposition using Feller's test for explosion. 

Let us start by considering a < 0. The diffusion is then defined on R with no singular point, 
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and Feller's scale function reads: 



p(x) = J exp |-2 £ -(? + a) l - 2a dfj dy. 



At infinity, the integrand inside the exponential is equivalent at infinity to 2(2Q , 2 1) _ 1 y~ 2(2g ~ 1)+1 . 
This quantity has a finite limit when a > 1, implying that — > ±00 when x — > +00. Feller's 
test implies that the first exit time of (-00, 00) is infinite with probability one. 

For a < 1, tends to 00 at 00, and has a finite limit at -00. This property ensures that 
the probability of the process diverge towards -00 is equal to one. It remains to prove that the 
explosion occurs in finite time. To this purpose, we need to show that the limit of Feller's scale 
function v(x) is finite at -00. This is the case under the assumption a > 1. Indeed, v(x) is defined 
as: 



r(.v) = £ £ CXV [ 2 .[ " 



The integral in the argument of the exponential function is equivalent at -00 to - z°) 

which tends to -00 as a power function, hence the exponential has moments of any order, and in 
particular is twice integrable at -00. 

For a > 0, we consider the intervals (-00, - y[d) and ( yfa, 00). The above analysis proved 
that for a > 1, p(x) tends to 00 when x — > 00 and towards a finite value at x — > -00. At 
+ yfa, the argument of the exponential term is equivalent to J x~ 2a+] which diverges for a > 1, 
implying that the exit time of (-00, - \fa), (- yfa, yfa) and ( \fa, 00) are all almost surely infinite 
and that there is no blow up in finite time. For a < 1, both p(- V« ) an d p( V^ + ) are finite. For 
c e (- y/a, y/a), we have: 

G(x) = ^) 1 - 2a (-^+ yfi?- 2a df 

, , . r r exp(-2(G(y) - G(z))) , , 

and using a similar argument as used in the proof of proposition 2j we show that lim t ^ ^ r v(x) < 
00, and by symmetry we hence have lim l .^_^+ v(x) < 00. Feller's test for explosion hence en- 
sures that the first exit time of the interval (- yfa, yfa) is almost surely finite. We hence converge 
in finite time towards yfa or - yfa. Moreover, the probability for a trajectory starting with an 
initial condition x e ( yfa) to reach y[a (resp. - y[a) prior to reaching - \[a (resp. \[d), denoted 
p+(x) (resp. p-(x)) have the expression: 

P-(x) = 1 - p+(x) - 



p{ Va ) - P(- Va + ) 

This quantity is plotted in Figure|9|c). □ 

Application of the general results of section [2] yield the following: 

Proposition 8. For a > 0, the saddle-node equation \\A) has two Dirac delta distributed func- 
tions 5 ^ and 6_ yg, which enjoy the stability properties: 

• for a < 1, both 6 ^ and S_ ^ are stable in probability, and the first hitting time t of one of 
the equilibria is almost surely finite. 
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• for a — I, the stationary distribution 5 ^ is always stable in probability. The stationary 
distribution S_ ^ is stable in probability for a > c~ 4 and unstable in probability otherwise. 

• for a > 1, the stationary distribution 6 ^ is stable in probability and is unstable in 
probability. 

Proof. This is a direct application of the results of section [2] applied to y + to characterize the 
properties of -\/a and to y~ for - yfa. For instance for the fixed point - -\/a and a = 1, the fixed 
point is stable if and only if 2 -\/a < 2a 2 a i.e. cr 2 y r a> 1. The finite absorption time property is 
a consequence of proposition [7] □ 

For a < 1, we know that the trajectories either reach one of the equilibria or blow up in finite 
time, which describe perfectly the permanent regime of the saddle-node bifurcation. We are now 
interested in stationary distributions of the saddle-node equation distinct from 6^ and 6_^. 
Since the system presents two singular points, the results of proposition|5] based on integrability 
properties on R + or Br of the solutions of the Kolmogorov equation, do not directly apply. 
However, a similar approach allows to demonstrate the following: 

Proposition 9. Stationary distributions of the saddle-node equation ( |14[ ) distinct from 6^ and 
(when a > 0) enjoy the following classification: 

• For a — 1 and < a < the system has a stationary distribution, stable in the first 
approximation, charging only the interval (— oo, — yfa) that undergoes a P-bifurcation at 
cr 2 y'a — j. On (— y/a, oo), any solution converges almost surely exponentially towards y'a 
(see Figure^. 

• For a — I and a < 0, there exists a unique distribution charging the whole real line. 

• For a > 1, the system has no additional stationary solution. 
Proof. In the case a = 1, the Kolmogorov equation has the solution: 

po(x) oc |x + Va| 2+lr2 ^ \x - V«| a? ~' ra 

For a > 0, this function is integrable at y'a if and only if -2-—^j= > -1 which is not possible. At 
- yfa, this function is integrable when -2 + > 1 i.e. cr 2 yfa < 1 (i.e. when - yfa is unstable). 
Eventually, this function is integrable at infinity whatever the parameters, since it behaves like 
x~ A . Hence there exists a unique stationary distribution charging (-oo, - y[a). In both the interval 
(- s/a, y'a) and ( yfa, oo) the only stationary distribution is 5 ^ which is stable. The additional 
stationary distribution on (-oo, - ^fa) reads: 

Po(x) = Z^-x - yfaj yfaj -:Zqo(x) 

Denoting <t> = ^r^, we give the following closed form primitives related: 



, w , r W+av r ,-i-<i>a(2<S> 2 - l)-2yfdd>x + x 2 

qo(x)dx = —(—ya — X) (ya-x) 



I 
I 

J~x 2 qo(x)dx = —{—y r a — x) l+ ®(~\fa-x) 1 ''' ^ 



4a3/ 2 <D(<D 2 - 1) 

, r \-i+av r w-o -2 rf*® x + a + xl 
x q (x) dx =(-ya-x) (ya- x) — — - — 

4a(<P z - 1) 



4 V^E>(<D 2 - 1) 
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Using these formulae, it is straightforward to obtain the normalization constant of p on (-00, - ya): 

Z = 4a 3/2 0(0 2 - 1) 

The expectation m is equal to: 



m ■■ 



xqoix) dx = -<t> a/a = ~ , 

o- 1 



and the second moment: 

X- V« 
x 2 q (x)dx = (2<t> 2 - l)a 

Eventually, let us analyze the shape of the obtained stationary distribution. The differential of the 
distribution reads: 



f J-2-9 -2 + 0\ 

P () (x) = po(x) + . 

\x — ya x + ya J 



This quantity is always negative for <£> < 2, the density is decreasing and moreover in that case it 
is clear that the distribution diverges at - yja. For <D > 2, the distribution vanishes at zero, takes 

its maximum at x* = -^5^ = - jo 3 m ^ 8 oes back to zero. 

Let us now address the stability in the first approximation for this distribution. The linearized 
equation reads: 

dv, = -2X,v, dt + 2o-X,v,dW t 

Integrating the linearized equation as previously done, we obtain that the Lyapunov exponent of 
the linearized equation reads: 

I = -2E[X t ] - 2cr 2 nxj] = 2{ ~ l + f a) 

(T 1 

and since we consider the case o 2 ya < 1, the Lyapunov exponent is strictly negative, imply- 
ing almost sure exponential stability of zero for the linearized equation, hence stochastic linear 
stability of po(x). 

For a < 0, the solutions of the Kolmogorov equation are given by: 

K / arctan(x/ V=a) 

(x l -a) 1 \ o 



and the definition interval of possible stationary distributions is R. In order for po to define a 
probability distribution, we need this function to be integrable at ±00, which is always the case. 
The obtained stationary distribution reaches its maximum again for x = - ■ 
For a > 1, the solution of the Kolmogorov equation reads: 

^=Ja^A^!i X{ - y2 + a) ^ dy ) 

At - yH, the exponential term diverges since it behaves as (x + y r a) 1 ~ 2a , hence the solution is 
never integrable at zero. It is easy at this point to show that there is no integrable solution to the 
Kolmogorov equation. □ 
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Remark. For a < 1 one could be interested in the existence of quasi- stationary solutions. An additional 
difficulty arises from the fact that the solutions can be attracted infinite time towards different solutions: 
± yfa whatever a < 1, and -co for a < |, Moreover, in the saddle-node bifurcation case as studied here, 
the diffusion coefficient is exactly equal to a power function. Therefore, for a < 1, blow up infinite time as 
well as absorption at ± ^Ja can occur are not a consequence of the results of U5V , are not in the scope of 
the present manuscript and require new mathematical developments that will be addressed elsewhere. 



(a) , (b) 




(e) (f) (d) 

Figure 9: Behaviors of the stochastic saddle-node system with a = 1. The bifurcation diagram as a function of a is at 
the center of the figure, and displays the singular points in blue lines, their stability marked by the type of curve: plain: 
stable, dashed: unstable. Stationary distributions are encoded in color in the interval a < <r~ 4 . Two typical stationary 
distributions are plotted in Figs, (e) and (f): one for a < l/4cr 4 where the pdf has a peak at x = -<r~ 2 /2, and for 
a > cr~ 4 /4 where the distribution diverges at zero. The P-bifurcation between these behaviors is noted P. (a) corresponds 
to a typical trajectory for a = 0.4 < cr~ 4 /4 with a = 0.8. (b),(c) and (d) correspond to a > cr -4 and different initial 
conditions: (c): frequency of the trajectories converging towards + ^[a as a function of the initial condition in [- yfa, ^fa\, 
(b) and (d) are typical trajectories for initial condition in ( V(a), oo) or (— oo, - yfa) respectively (a = 4 and cr = 0.8). 
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The above result characterize the dynamics relatively exhaustively for a = 1, except for 
a > cr~ A and xq € (- ^[a, yfa). In that case, both singular points are stable, and can be reached 
by the solution. Therefore for such initial conditions, there is a competition between the possible 
attractors, and the relative stability of one singular point compared to the other will determine 
the probability, given an initial condition xq, to converge towards this equilibrium. This quantity 
is plotted in Figure j^c). We observe, as one can expect, that the attractivity of Va is stronger 
than that of - -\/a: the probability of converging towards -\/a is never less than 0.4 even for initial 
conditions very close from - -\/a, and goes relatively fast towards 1 as the initial condition is 
chosen closer from yfa. 

Let us now discuss the universality of this bifurcation. We consider a general one-dimensional 
diffusion equation: 

dx, — f(x,, A) dt + g(x,)dW t 

with / a smooth function such that /(0,0) = §j(0,0) = and §^(0,0) + 0. Following the 
reduction to a normal form of Kuznetsov (2), we use a Taylor expansion of /: 

f(x, A) = f (A) + fx (A) x + MX) x 2 + Kx) 

where h(x)/x 3 is bounded at x = 0. We perform a shift of coordinates by defining a new variable 
£ = x + 6 and obtain: 

= {[MA) - Mm + fi{X)6 2 + 0(6 3 )] 
+ [f 1 (A)-2f 2 (A)6 + 0(6 2 )^ 
+ [f 2 (A) + 0(6)] f + 0{f)) dt + - S)dW, 

Using the inverse function theorem under the assumption that |j(0,0) + 0, we deduce that there 

exists a smooth function 6(A) = ^q^A+0(A z ) such that the term [/[ (/l)-2/ 2 (/l)(J+(9((5 )] vanishes 
for all A in a neighborhood of 0. This 6(A) identified, we can apply the inverse function theorem 
to find locally an inverse to the constant term of the Taylor development seen as a function of A. 
This inverse function, A(jS), yields the system to the SDE: 

d& = Ou + b(^] + 0(f)) dt + g(g - 5(A(ji)))dW, 

with b(Q) = M®) ^ 0- The final reduction involves a scaling of the variable: z, = \b(p)\^ t , 
yielding the SDE: 

dzt = (p + sz] + 0(z 3 )) dt + g(-^ - 6(A( M )))dW, 

where e has the sign of 5^(0). The stability and bifurcations of the whole system is hence 
reduced to the local analysis of the behavior of g close of the equilibria. 

Example 3. We illustrate this method using a specific instance of a one-population Wilson and 
Cowan system with constant external input I e R: 

dv, = ( - v, + S (gv t ) + l)dt + o-\-v,+S (gv,) + I\ a dW, 

and choose as usual S (x) — tanh(g x). The whole dynamics is hence governed by a single function 
f(x) = -x + S(gx) + I and the deterministic system undergoes two saddle-node bifurcation when 
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v — +— — 



g > 1, located at v — + — /or an /n/?Mf equal to ( respectively). 



I = ±^E±+S(J^T). Let 



us focus for instance on the bifurcation arising at 
yields, denoting x t -v,- 



. Around this point, a Taylor expansion 



dx, = 



/- 



+ 5(V^T) - J^Txj+h(x t ) 



dt 



+ o- 



g 



1 ) 



4g^\x 2 t + h(x,) 



dW, 



where h(x) = 0(x 3 ) at zero. Eventually, denoting jj t = Jg — 1 jc r we obtain: 

= g + dt+\p-g+ vm a dw t 



with p - Jg — 1(7 + 5 ( - 1) — - — anrf = -v/g - 1 /z( p — ). 77ie system is hence locally 

reduced to the saddle-node equation, and the above analysis characterizes the dynamics of the 

system. Let us for instance consider I as our parameter. As I crosses I* — ~— — S ( Jg — 1), 
a stochastic saddle-node bifurcation occurs. For I > I* the system has two equilibria xi < 
Jg — 1 < x r and an additional singular point x* strictly smaller than x/. It is easy to show that 
this singular point is always stable in probability. Regarding xi and x r , we have: 

• For a < 1, both x% and x r are stable and the first hitting time of one of the equilibria is 
finite almost surely. 

• For a = 1, x r is always stable in probability, and x/ is stable for (3 > o~~ A , and unstable 
otherwise. 

• For a > 1, x r is stable in probability and xi unstable. 



We illustrate this variety of behaviors in Figure 10 



4. Discussion 

In this manuscript we analyzed the interplay between noise and nonlinearities in the setting of 
SDEs with singular points, and showed that the present analysis could account for the dynamics 
of a number of low-dimensional stochastic systems driven by a unidimensional Brownian motion. 

Even though these limitations restrict the variety of behaviors possible, we observed and 
quantified several effects of noise, such as for instance the that it can prevent explosion of the 
solutions, and that it interferes with the stability of singular points. Moreover, as soon as the 
dimension of the vector field is increased (e.g. Hopf bifurcation), stochastic perturbations can 
have more complex effects in which the precise structure of the noise interacts subtly with the dy- 
namics. In particular we identified a new bifurcation corresponding to a change in the orientation 
rotation of the stochastic cycles. 

An important question related to this manuscript is the possibility to reduce the dynamics 
of general SDEs to low-dimensional systems, analogous to the center manifold reduction of 
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(a) a = 0.8, <x = 2 (b) a = 0.8, <x = 2 (other sample (c) a = 1, cr = 1 

path) 




(d) a = 1, <r = 2 (other sample (e) a = 1, cr = 0.2 (f) a- = 1.5, cr = 0.2 

path) 

Figure 10: Dynamics of the Wilson and Cowan System (v, as a function of time) for different values of a in the saddle- 

node case with g = 3, / = 0, and initial condition - s — — . In that case, x r = and x\ has no closed-form expression. We 
illustrate the stabilization of xi in the case or < 1, which holds only for large noise values when a = 1: (a) and (b): a < 1, 
.r/ and a> are stable in probability and the first hitting time to one of the equilibria is a.s. finite, (c) a = 1, cr large enough 
so that both equilibria are stable, (d): for cr small, only x\ is stable, (d): for a > 1, only x\ is stable. The initial condition 
is taken in the interval (x r , X(), and the stability of .v* is not addressed. 
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deterministic dynamical systems. This is highly dependent on the type of diffusion coefficients 
considered, and the dimension of the Brownian motion involved in the dynamics. Indeed, in the 
stochastic case, when performing nonlinear changes of variables, additional terms appear in the 
drift function, that depend on both the drift and diffusion coefficients. This analysis is part of 
the perspectives of this work. Another perspective is the further analysis of the quasi-stationary 
distributions of such systems. Indeed, the present analysis demonstrated that quasi-stationary 
distributions existed as long as a < |. The question hence arises whether this finding is a 
limitation of the framework we used, or if actually disappears and we are left with no stationary 
behavior for the trajectories that did not reach zero. This would be the first kind of process 
discovered, up to our knowledge, of diffusion absorbed in finite time, with no quasi-stationary 
distribution. We believe that in the case A > 0, the quasi-stationary distributions progressively 
turn into stationary distributions as the exponent a increases to 1. This would correspond to 
a convergence towards zero of the eigenvalue of the Kolmogorov operator related to the quasi- 
stationary distribution. But for A < and a > 1, there is no stationary distribution for the process, 
and hence these quasi-stationary distribution do disappear. And numerical simulations provided 
in Figure [2] tend to confirm that indeed, the quasi-stationary distribution disappear in favor of 
a progressive convergence towards zero of all trajectories conditioned on not hitting zero. The 
analysis of this phenomenon is another perspective for future work. 
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